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Parametric  Polynomial  Curves  of 
Local  Approximation  of  Order  8 


K.  Scherer 


Abstract.  Parametric  approximation  of  curves  offers  the  possibility  of 
increasing  the  order  of  approximation  by  using  the  additional  parameters 
in  the  parametrization  of  the  curve.  This  has  been  studied  in  several 
papers,  see  e.g.[l-8].  The  resulting  problems  are  highly  nonlinear.  Here 
the  cases  of  approximation  order  0(h&)  are  studied  which  need  piecewise 
quartic  curves  in  the  plane  and  piecewise  quintic  curves  in  space. 


§1.  Introduction 

A general  conjecture  concerning  the  local  approximation  order  by  polynomial 
curves  can  be  formulated  as  follows  (see  e.g.  Mprken-Scherer  [6]): 

• For  sufficiently  small  h > 0 and  a sufficiently  smooth  curve  f(t)  : t £ 
[0,  h]  — > /(t)  € Rd,  there  exists  a polynomial  curve  p(t)  of  degree  n and 
a reparametrization  ip  of  / on  [a,  6]  such  that 

sup  ||(/  o <p)(t)  - p(t)||  < C(f)hm  , 

0 <t<h  ~ ~~ 

where  m :=  n + 1 + [*j5y]- 

The  increased  order  m is  explained  in  [6]  by  the  principle  of  degree  reducing. 
It  comes  from  the  idea  of  approximating  with  an  interpolating  polynomial 
curve  p(t)  of  degree  m — 1 such  that 

p{U)  = (/°  4>){ti),  1 <i<m,  (1) 

for  points  f,  in  [0,  h]  (multiplicities  allowed).  The  additional  parameters  oc- 
curing  via  the  reparametrization  <j>  are  used  to  reduce  the  degree  m—  1 of  p(t) 
by  requiring 

[ti,...,tm_i](p)  = [fi,...,fm_j](/o</>)  = 0,  « = 0,  1,  (2) 
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with  k as  large  as  possible.  Since  we  have  to  normalize  0 such  that  0(0)  = 
0, 4>{h)  = h,  there  are  m — 2 parameters  left  at  our  disposal  for  this  goal.  Thus 
m— 1 can  be  reduced  ton  = m—l  — k,  where  k-d  < m — 2 or  k < (n  — \)/{d—  1). 
Prom  classical  approximation  theory  it  is  then  clear  that  a solution  of 

(1) ,(2)  yields  a polynomial  curve  p(t)  of  degree  m—l—k  satisfying 

sup  \\(f  o - p(t)\\  < hm  sup  \\Drn(f  o 0)(t)||.  (3) 

0 <t<h  0 <t<h  ~ 

Thus  conjecture  (•)  is  true  if  one  can  guarantee  in  addition  that  the  parame- 
ters 0(2)(O), . . . ,0<m)(O)  of  a solution  of  (2)  remain  bounded  for  h — * 0.  This 
question  of  stability  is  also  discussed  in  [6].  Note  also  that  in  this  case  equa- 
tions (1)  can  be  written  as 

pO0_1(sj)  = /(fii),  1 < 7 < 771, 

where  the  nodes  Si  are  defined  by  s,  :=  0(t,). 

The  most  interesting  case  of  the  conjecture  is  when  (71  — 1 )/{d  — 1)  is  an 
integer  k,  i.e.  77i  = n + l-l-A:  = fc<i  + 2.  In  this  case  there  are  kd  equations  in 

(2) ,  and  the  degree  771  — 1 is  reduced  by  k.  Then  0 is  determined  by 

[ti,...,t„+i+i](/o0)  = 0,  i = (4) 

and  p(t)  by  the  first  n + 1 equations  in  (1). 

So  far  the  conjecture  seems  to  be  proved  only  in  the  case  n = d or  k = 1 
(see  [5-6])  and  for  k — 2 = d (see  [1,3]).  Here  we  treat  the  next  most  difficult 
cases  k = 3,d  = 2 and  k = 2,d  = 3,  which  amount  to  six  equations  in  (4), 
and  will  lead  to  quartic  curves  in  the  plane  and  quintic  curves  in  space  with 
approximation  rates  of  order  8,  respectively. 

§2.  Reduction  to  2 x 2 Systems 

In  Morken  -Scherer  [6]  equations  (4)  were  studied  in  particular  for  the  Taylor 
case 

Dn+i{f_  o 0)  = 0,  t = l,...,Jb. 

The  crucial  point  is  then  the  formula  of  Faa  di  Bruno  in  the  form  given  by 
T.  Goodman.  It  reads  (cf.  [4]) 

( 

Dl(loct>)(0)  = Y,ai,dU](Q),  Pi  ■=  0(O(O),  (5) 

i= 1 

where 

ai,j  = y~i 

h+...ij=i,i>  1 

1 1 n 

h”‘lj\  Ip.  ■ ■ ■ Ijlmp.  ■ ■ ■ mr\ 


h--lj 


Ph'-Plj 


and 
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Here  the  integer  r denotes  the  number  of  distinct  integers  among  the  li, . . . , lj, 
where  mi,...  ,mr  are  the  multiplicities  of  them.  Specific  examples  are 


ai,i  = Pi, 


«i,i  = Pi- 


At  first  we  consider  the  planar  case  d — 2 for  k = 3.  Then  the  equations  (4) 
specialize  to 

Di(foct>)(0)  = 0,  j = 5,6,7.  (6) 

Under  the  normalization  <j>( 0)  = 0,  Pi  :=  <f>'{ 0)  = 1,  this  yields  six  equations 
for  the  unknowns  P2,...,Pt-  In  [8]  this  system  has  been  reduced  to  a 2 x 2 
system  for  P2,P^-  The  basic  idea  was  to  simplify  system  (6)  by  determining  a 
preliminary  reparametrization  ip  with  ip( 0)  = 0,  ip'( 0)  = 1 such  that 


•D2i+1(/°V 0(0)  =0,  i — 1,2,3. 


It  can  easily  be  shown  that  this  system  is  uniquely  solvable  in  the  unknowns 
72,  ...,77,  where  7,-  :=  0)  provided 

span(/'(0),/"(0))  = lR2.  (7) 

Thus  one  can  assume  that  D3f( 0)  = D5/(0)  = D7f( 0)  in  (6).  Denoting  dij 
as  the  cross  product  of  fU)  and  /W  in  1R2,  i.e. 

dij  :=  f}»  x /(i). 

Further  straightforward  computation  (cf.  [9]  ) leads  to 

Lemma  1.  Under  the  assumption  (7)  and  the  normalization  (j>( 0)  = 0,  Pi  := 
<j>'(0)  = 1,  the  vector  P2,  ■■■,p7  is  a solution  of  the  equations  (7)  iff  P2,  P3  is  a 
solution  of  the  2x2  system 


0 — ^6,1  + 20(^47^3  + 10^2,1^3 

+ 60p2d4,2  + 15(^4, 1 — 2d2,\Pz)P2, 


(8) 


and 


0 = 15/3^4,1  + 75/32d4,2  + (3dg,i  — 10d%1/d2li)P2  + ^6,2  . . 

+ 10(2d4,2  + P2d4,i)P3  — Wd2,iP2Pl 

The  elimination  of  the  variable  ps  in  this  system  along  the  lines  of  the 
resultant  method  yields  an  equation  of  degree  9 in  /?2-  However,  it  was  over- 
looked in  [8]  that  the  constant  term  in  this  equation  vanishes,  so  that  it  is  in 
essence  of  degree  8 . Thus,  existence  of  a (real)  solution  of  (8)-(9)  cannot  be 
derived  in  this  way.  We  will  close  this  gap  in  the  next  section. 
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In  the  space  case  d = 3 we  have  to  find  a reparametrization  <f>  such  that 

Dj(£o<t>)(  0)=0,  j = 6,7.  (10) 

We  simplify  this  again  by  determining  a reparametrization  with  i/j(0)  = 0, 
ip'(0 ) = 1,  ip^(0)  :=  7 j and 


Dj(£  oV-)(0)  = 0,  j = 4,7. 

This  is  possible  since  by  (5),  in  the  three  equations  forming  Di(f  o ip)( 0)  = 0 
the  coefficients  of  f',f"  and  f”'  are  linear  in  the  unknowns  72,73,74,  and  in 
the  latter  three  are  linear  with  respect  to  the  75,  76,  77.  Therefore,  these 
equations  are  uniquely  solvable  under  the  assumption 

span  (/'(O), /"(0), /"'(0))  = IR3,  (11) 

and  we  can  consider  (10)  without  loss  under  the  assumption  D4(/)(0)  = 
D7 (f)(0)  = 0.  This  yields  (with  0i  = 1)  the  equations 

0 = 06/'  + (605  + 150204  + 10/?|)/" 

+ (15  04  + 6O0203  + 15/?|)/"'  + 1502/(5)  + /(6),  1 

0 = fhf  + (70s  + 21020s  + 350304)/"  + (1O5022  + 35  03)/(5) 

+ (2105  + 1050204  + 7001  + 1O50203)/w  + 2102/(6).  ' ^ ' 

The  next  step  is  to  take  in  (12)  the  scalar  product  with  cross  products  }'  x 
f"  x f"'  and  f"  x f"' , respectively.  We  obtain  the  equivalent  equations 

0 = (1504  + 600203  + 1501)^1,2,3  + 1502^5,1,2  + <^6,1,2  (14) 

0 = (605  + 150204  + 1001)^2,1,3  + 1502^5,1,3  + ^6,1,3  (15) 

0 = 06^1,2,3  + 1502^5,2,3  + ^6,2,3,  (16) 

where 

(/(il)  X /(<2),/(j))  :=  det(/0'),/02),/0))  ~ diuhJ. 

These  equations  serve  for  eliminating  the  variables  04,05  and  06  since  they 
appear  linearly.  Before  doing  this,  we  transform  (13)  into  three  equivalent 
scalar  equations  analogously  to  (12).  We  obtain  the  three  equations 

0 = 07^1,2,3  + (10501  + 3503)^5,2,3  + 2102^6,2,3,  (17) 

o = (06  + 30205  + 50304)^2,1,3  + (1502  + 503)^5,1,3  + 302^6,1,3,  (18) 

0 = (305 + 150204  + 1003 + 150203)d3,l,2  + (1502 +503)^5, 1,2  + 302^6, 1,2-  (19) 
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Equation  (17)  determines  4 directly  in  terms  of  4 >4-  Now  we  eliminate 
4>4  in  (18)  via  (15), (16).  This  gives  at  first 

o = — (344  + 544)4,1,3  + (15/32  + 54)4,1,3  + 34(4,1,3  + 54,2,3)  + 4,2,3 

and  then 

15  5 

0 = 02^1,2,3  + 544,1,3  + 2 ^2(^6, 1,3  + 64,2,3)  + 4,2,3 

15 

+ (~2~/324  + 5/3f4)4,2,3  — 5444,2,3- 

Now  we  eliminate  the  variable  @4  by  (14).  The  result  is  the  equation 

0 = 9o(/32)  + 4(4)4  + 92  (/32  )/?3 » (20) 

where 

4(4)  :=  (15/2)4, 2, 3/325  - ( 15/2)4, 2, s/323  + [(15/2)4,3,5  - (l/2)4,2,6]/322 

+ (154,3,5  _ (5/2)4, 3, 6)4  + 4,3,6) 

4(4)  :=  (1/3)4, 2, 6 + 54,3,5  + 54,2,54  — 254,2,3/32> 

4(4)  :=  —254,2,34- 

Analogously  we  reduce  equation  (19)  to  an  equation  in  4>4  by  eliminating 
4 and  then  4-  Using  (15),  we  obtain 

0 = (5/3|  + 15/3|4)4,2,3  + (15/3|  + 54)4,1,2 

+ 34(4,1,2  + 24,1,3)  + '~6'2  '3  3"  "2’444,2,3> 

and  then  by  (14) 

0 = 4,2,3(— 15/3|  + 30442  + 1042)  + 4,1,2(15$  + 104)  ,21) 

+ 4(154,1,3  + 54,1,2)  + 4,1,3- 

In  order  to  get  rid  of  the  term  with  /3|  in  (21),  we  make  a final  substitution 

4 :=  4 + a/3|,  a :=  (3/2)  + y/IB/2. 


Then 

-1544  - 30442  + 1042  = [-3042  + 20a/3f]4  + 1042, 

and  (21)  simplifies  to 


0 = po(4)  + pi  (4)4  + P2(4)/3|, 


(22) 
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where 

P0W3)  10^1, 2, 3ft2  + 10di,2,5^3  + ^1,3,6, 

Pi(/%)  :=  5di,2,6  + 15di,3l5>  (23) 

P2(ft)  :=  10\/i5rfl,2,3^3  + (15  + 10a)di,2,5- 
With  the  new  variable  P3,  (20)  transforms  into 

0 = 90(^2)  + 9\{P2)Pz  + 92(^2) 032 , (24) 

with 

Qo (/?2 ) :=  9o(&)  - 25a2rfi,2,3,3|  + a/3 lqi(P2), 

9l(^2)  :=  9l(ftf)  — 50^1,2,30^2) 

92(^2)  :=  92 (P2)  = — 25di,2,3/?2, 
and  the  qi(P2)  defined  as  above.  We  summarize  all  this  in 

Lemma  2.  Under  the  assumption  di,2,3  i1  0,  i.e.  assumption  (11),  and  the 
normalization  0(0)  = 0,/3i  :=  0'(O)  = 1,  the  vector  p2,  is  a solution  of 
the  equations  (22),  (24)  iff  p2,  (h  is  a solution  of  the  2x2  system  (22),  (24). 

Remark:  The  systems  (8)-(9)  in  the  planar  case,  and  (22)-(24)  in  the  space 
case  possess  a similar  structure.  In  (8)  and  (22)  the  coefficients  of  P2  are 
polynomials  of  the  same  degree  in  p 3 and  p2,  respectively.  The  same  is  true  for 
(9)  and  (24),  except  that  the  corresponding  polymials  have  different  degrees. 


§3.  Existence  Theorems 

In  view  of  the  last  remark,  we  treat  in  detail  only  the  planar  case. 

Theorem  1.  The  system  (8)~(9)  has  at  least  one  and  at  most  5 (real)  solution 
pairs  P2,Pz  outside  the  line  p2  = d4,i/2d2,i  ■ 

Proof:  Let  us  write  for  shortness  x :=  P2  and  y :=  P2  as  well  as 


A(y)  :=  Po(y),  2 B :=  pi(y)  = 60 <f4,2,  C{y)  :=  p2(y)  = 15d4,i  - 30d2,iP- 
Then  (8)  reads  0 = A(y)  + 2 Bx  + C(y)x2.  Formal  solution  for  x gives 


= <p±(y)  := 


-B  ± sjB2  - A(y)C(y)  _ - B ± y/R&j 


C(y) 


C(y) 


(25) 


with  the  cubic  polynomial 


R(y)  — 15[20ci|ji  j/3  + 30d4]id2,ii/2  + (2(16,1^4,1  — 20d4il)i/  + 240d22  — 1 


6, 1 J 


Then  write  (9)  as  0 = Y^=o  ttiX'  4"  + b\xy  + b2xy2,  and  insert  (25),  since 

by  assumption  C(y)  ^ 0.  After  multiplication  with  C(j/)3,  we  obtain 


0 = X>(-B  ± + b0yC(y)3 


i=o  (26) 

+ h(-B  ± y/R(y))C(y)2  + b2y2(-B  ± J W))C(y )2 ■ 
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Now  observe  that 

(-B  ± v^))2  = B2  + R(y)  ± 25^^), 

(-5  ± \/«(y))3  = - B 3 ± W2y/R(i)  - 3BR(y)  ± R(y)y/R^), 

and  sort  all  terms  with  and  without  \/R(y),  respectively.  Then  (26)  can  be 
written  as 

U(y)  = ±V(y)y/RG),  (27) 

where  U(y)  is  a polynomial  with  leading  term  —11  • 15  • ■ 9000  d4^  d\  \ y4  and 
V(y)  also  of  degree  4 with  leading  term  9000  d%  4 y4. 

Hence  under  the  above  assumption,  P2,P3  is  a solution  of  the  system 
(8)-(9)  iff  y = /%  is  a solution  of  (27)  with  sign  either  + or  - on  the  right 
hand  side.  Suppose  that  n\  and  ri2  are  the  numbers  of  solutions  of  these  two 
equations  (including  multiplicities).  Then  the  squared  equation 

U\y)  = V\y).R{y) 

is  of  degree  11,  and  has  2(ni  + 712)  solutions.  Hence  we  conclude  n\+U2  < 5. 
To  prove  existence,  write  (9)  as 

H(/32,p3)  = 0,  (28) 

where  H is  of  degree  3 in  [i2  and  with  -10/32/?|d2,i  as  leading  term  in  f33.  Then 
introduce  y*  as  the  largest  zero  of  R(y),  so  that  in  view  of  R(+ 00)  = +00 

R(y*)  = 0,  R{y)  > 0 for y*  < y < 00. 

Now  insert  both  functions  ip±(y)  in  (28),  and  obtain  the  functions 

H+{y)  :=  H(ip+{y),y). 

In  order  to  guarantee  existence  of  a solution  of  (8)-(9),  it  suffices  therefore  to 
show  that  the  ranges  of  H+  satisfy 

H+[y*,  00)  U [3/*,  00)  = IR.  (29) 

For  this,  observe  at  first  the  properties 


My*)  = -B/C(y*)  = 2d4,2/(d4,i  - 2 y*  d2,i)  :=  <P* 


and 

P+(y)  « ±i||1/2sign  (d2,i),  y-KX). 
Then  we  distinguish  the  cases  (assume  without  loss  d2,i  > 0): 


i)  y*  > y :=  diA/2d2ti, 


ii)  y*  < y. 
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In  case  i),  we  consider  the  ranges  of  ip+  for  (y* , oo),  and  have 

V+[y*,°°)  = (-0O,  v>*],  ¥>-[y*,oo)  = [<p*,oo).  (30) 

Since  both  tp+  are  well  defined  and  continuous  on  (?/*,  00),  so  are  the  functions 
H+(y),  and  furthermore 

H±{y*)  = H(<p*ty*). 

In  combination  with 

H±(y) « 10^2, i||3/|2v+(y).  y ->  (31) 

it  follows  that 


H+[y*,oo)  = (—oo,H((p*,y*)},  H.[y*t  00)  = [#(?*,»*),  00), 

and  hence  the  desired  assertion  (29)  in  case  i). 

In  case  ii),  the  function  <p~(y)  has  a singularity  at  y which  lies  in  (y* , 00). 
However,  we  can  restrict  its  domain  to  [y*,y)  and  still  have  (d2, 1 > 0) 

•p-[y*>y)  = [<?*>°°)-  (32) 

On  the  other  hand  the  function  ip+  (y)  remains  continuous  on  the  whole  inter- 
val (y*,  00)  since 

<P+(v)  = MV)/{B  + VB2  - A(y)C(y)). 


Hence  we  have 

ip+[y*,oo)  = (—00,  v?*], 

so  that  together  with  (32)  we  have  the  same  situation  as  in  (30)  and  can 
proceed  further  exactly  as  before  in  order  to  prove  (29).  □ 

It  remains  to  discuss  whether  there  exist  solutions  of  (8)  (9)  if  /33  = 
di, 1/2^2, 1 :=  y ■ In  this  case,  (8)  gives  /02  = A(y)/B  if  B / 0,  and  (9) 
can  have  a solution  only  under  some  additional  constraint  on  the  parameters 
d2,ii c?4,i , ^4,2, ^6,1) ^6,2-  We  omit  it  here,  as  well  as  the  one  which  results  from 
(8)  if  in  addition  B = 60d4i2  = 0. 

Further,  we  remark  that  there  can  indeed  exist  5 solutions  of  (8)-(9).  To 
this  end,  one  can  consider  the  case  dji  = 0,  where  these  equations  simplify 
in  such  a way  that  solving  (8)  for  z :=  /? 3 — (3/2)/3|  gives 


±z 


2^4, 2P2  1 
^2,1 


de,ip2  2 

30d2  ! 


+ 0(/?2~4). 


Inserting  this  into  (9)  with  sign  +,  it  follows  that 


0 — — 90d2,i/^2  T 255(14,2(02  T 5d6,i/?2  + ^6,2  + 0(1 /fa)- 


Here  the  polynomial  of  degree  5 dominates  for  large  /?2,  and  it  is  clear  that  the 
4 parameters  <^4,2, o?6,2, <^6,1  can  be  chosen  such  that  there  exist  5 zeros 
outside  some  bounded  interval  containing  0.  The  situation  for  the  system 
(22)-(24)  in  the  space  case  is  similar. 
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Theorem  2.  The  system  (22)-(24)  has  at  least  one  and  at  most  7 (real) 
solution  pairs  ft,  ft  if  ft  / -(6  + v^)ft,2,5/5\/l5ft,2,3- 

Proof:  Concerning  existence,  the  argument  is  the  same  as  in  Theorem  1.  We 
define  x :=  ft,y  :=  ft  and  A(y)  :=  p0(y),B  :=  pi{y),C(y)  :=  p2(y).  Then 
R(y)  in  (25)  has  the  form 

R(y)  = 50ftl5  d\xzyz  + 50ftl5ft,  2,5 y2  + linear  term. 

Again  define  y :=  —(6  + v^)ft,2,5/5ftl5ft,2,3  as  the  zero  of  C(y)  and  y*  as 
the  largest  zero  of  Aft),  and  let  <p±(y)  :=  (B  ± ftAft))/Cft).  Its  asymptotic 
behaviour  is  described  by 

P±(y)  «sgn  dit2,3\-^\1/2<  y oo.  (33) 

Now  we  distinguish  as  in  Theorem  1 the  cases  i)  and  ii)  and  conclude  that 
either  (30)  holds  or  (31),  respectively. 

Next  write  (24)  similarly  as  in  (28)  as 

£ (ft,  ft)  =0 

and  define  fttft)  :=  H(tp±(y),y).  Its  asymptotic  behaviour  is  somewhat 
more  complicated  to  determine  than  in  (31),  since  by  the  definition  of  the 
<7ift),  2 = 0, 1,2,  in  (24),  the  leading  term  of  ii  (ft,  ft)  now  has  the  form 

( ~2 — 25a2)ft,2,3/3f  + 25aft,2,3/32  + (25ft,  2, 3 — 50aft,2,3)ft/3f  — 25ft,2,3ftft2. 
Inserting  (33)  with  ft  — y,  one  derives  from  this 

H±(y)  « const.  \y\5^2,  y — > 00. 

Hence  we  obtain  the  same  property  (29)  for  H±  as  for  H±  in  Theorem  1,  and 
existence  of  a (real)  solution  pair  of  (22),  (24)  follows. 

In  order  to  show  the  bound  on  the  number  of  solutions,  we  proceed  again 
as  in  Theorem  1.  We  solve  (22)  for  x = ft  as  in  (25)  and  insert  it  into  (24) 
written  as 

5 3 

0 = ajx'  + y y~]  bjx'  + cxy2, 

i=0  i=0 

with  constants  ai,bi  and  c.  This  gives  after  multiplication  with  C(y)5 
5 

0 = '£<p(-B±Vm)iC(y)5-i 

i=0 

3 

+ V h{-B  ± x/^))iC(y)5-<  + cy\-B  ± yfR^))C(y)\ 

i=0 

Then  sort  by  terms  with  and  without  ± ft  Aft)  . The  result  is  an  equation 
of  type  (27),  this  time  with  polynomials  U(y)  of  degree  7 and  V(y)  of  degree 
6.  Therefore  we  conclude  by  the  same  argument  as  in  Theorem  1 that  the 
system  (22),  (24)  has  at  most  7 solutions  . □ 

The  discussion  of  the  degenerate  case  ft  ^ — (6  + ftl5)ft,2,5/5ftl5ft,2,3 
is  omitted.  It  can  be  done  similarly  as  for  Theorem  1. 
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§4.  Final  remarks 

We  have  shown  that  local  approximation  order  8 can  be  achieved  with  para- 
metric polynomial  curves  of  degree  4 in  the  planar  case,  and  degree  5 in  the 
space  case.  The  method  for  this  consists  in  determining  a suitable  reparame- 
trization  <f>,  and  then  the  Taylor-polynomial  with  respect  to  fo<j>.  For  practical 
purposes,  however,  it  is  important  to  consider  also  Lagrange  or  Hermite  inter- 
polation in  the  sense  of  the  equations  (1).  This  has  been  done  in  case  k — 2 
for  d = 2 in  [1,3,6]  and  for  d = 3,  in  [5],  but  results  for  higher  k do  not  seem  to 
be  available  so  far.  In  this  respect  another  interesting  open  question  is  which 
order  of  geometric  continuity  can  be  preserved  when  a piecewise  polynomial 
curve  is  constructed  by  pieces  of  such  local  aproximations. 
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